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Abstract 

The tensor model can be regarded as theory of dynamical fuzzy spaces, and 
gives a way to formulate gravity on fuzzy spaces. It has recently been shown that 
the low-lying fluctuations around the Gaussian background solutions in the tensor 
model agree correctly with the metric fluctuations on the flat spaces with general 
dimensions in the general relativity. This suggests that the local gauge symmetry 
(the symmetry of local translations) is also emergent around these solutions. To 
systematically study this possibility, I apply the BRS gauge fixing procedure to 
the tensor model. The ghost kinetic term is numerically analyzed, and it has been 
found that there exist some massless trajectories of ghost modes, which are clearly 
separated from the other higher ghost modes. Comparing with the corresponding 
BRS gauge fixing in the general relativity, these ghost modes forming the massless 
trajectories in the tensor model are shown to be identical to the reparametrization 
ghosts in the general relativity. 
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1 Introduction 



Various thought experiments considering quantum gravitational fluctuations have shown that 
the classical concept of smooth spacetime in the general relativity is not appropriate in some 
extreme cases [TJ [2], and should be replaced in some way by a novel concept of quantum 
spacetime. Fuzzy spaced is one of such candidates of quantum space [31 El E]- A fuzzy space 
is defined by an algebra of functions on it, unlike a classical spacetime being described by a 
coordinate system. This kind of algebraic definition of spaces has some physically interesting 
advantages over the classical description. For example, in quantum gravity, the changes of 
topologies and dimensions of space are believed to be the vital processes of quantum fluctua- 
tions. However, it is generally hard or tightly constrained to describe these processes without 
encountering singularities in the classical description [7J. On the contrary, in general, one 
will have much more freedom to describe such processes in fuzzy spaces through interpola- 
tion between algebraic structures of fuzzy spaces approximating classical spaces with distinct 
topologies and/or dimensions. This kind of thoughts suggest an interesting research direction; 
considering theory of dynamical fuzzy spaces as a model of quantum gravity. 

In the recent years, there have been numerous discussions about gravity on fuzzy spaces. 
A class of approaches discuss analogues of the general relativity on fuzzy spaces. In this 
class of approaches, however, the dynamical variable is a fuzzy analogue of the metric tensor, 
and a fuzzy space itself is assumed to be fixed. Therefore these approaches do not take the 
full advantages of the notion of fuzzy space as explained above. On the contrary, a more 
interesting kind of approaches were initiated by the matrix models [Bl E]. These approaches 
consider spaces as dynamical objects generated as classical solutions or vacua, and fluctuations 
of matrices around such vacua are regarded as field fluctuations on background fuzzy spaces. 
Then an extremely interesting possibility is that gravity may appear as one of these emergent 
fields. So far this is yet an open issue under active investigations [101 [11] . 

In view of this present status, it might be meaningful to study another kind of model of 
dynamical fuzzy spaces, which is similar to but distinct from the matrix models. The model 
studied in this paper is the tensor model, which has a rank-three tensor as its dynamical 
variable, instead of matrices in the matrix models. The tensor model was originally proposed 
to describe the simplicial quantum gravity in dimensions greater than two ^ E31 M 113 ESI EH 
118)1 . The tensor model has not yet been successful for the analysis of the simplicial quantum 




*In this paper, this term is used as its widest meanings. It includes noncommutative spaces as well as 
nonassociative ones [3] . 

tightly related kind of models, called the group field theories, have been being discussed mainly in 
the context of the loop quantum gravity. It is known that a certain group field theory can be considered 
to be a field theory on a noncommutative spacetime [19] and can also be derived as effective field theory of 
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gravity itself, in part because of the absence of the analytical methods to solve the tensor 
model. However, it was recently proposed by the present author that the tensor model may 
be reinterpreted as theory of dynamical fuzzy spaces [231 1211 125]. This is based on the fact that 
a fuzzy space can be characterized by a rank-three tensor C a ^ c which determines the algebraic 
relations among all the functions f a on a fuzzy space through the product f a * fb = C a b c f c - 
From this point of view, it may not be necessary to analytically solve the tensor model to make 
relations to physics. In analogy with the matrix model mentioned above, a classical solution 
in the tensor model may be regarded as a background fuzzy space, and the fluctuations of 
the tensor about it as field fluctuations. Then the question is whether gravity appears in 
such fluctuations. In fact, in a class of tensor models which have the classical solutions with 
Gaussian forms, it has been shown that the low-lying fluctuations about these solutions at low 
momenta match correctly with the metric fluctuations on flat spaces in the general relativity 
in general dimensions (26J, [271 EE] • 

The above agreement is very interesting, but it is merely classical and obviously not enough 
quantum mechanically. The main purpose of this paper is to show the agreement a step further 
to include the gauge degrees of freedom, which are the local translations in the present case. 
The tensor model has the symmetry of the orthogonal group O(N), where N is the number 
of all the functions or more physically "points" forming a fuzzy space. A background solution 
of the tensor model breaks this O(N) symmetry down to some remaining symmetries of the 
solution or the background space, and the broken symmetries are realized non-linearly around 
hj*. Since the broken symmetries permute the "points" of the background fuzzy space, they are 
intrinsically local symmetries, and it is tempting to insist that these are emergent local gauge 
syrettes (the local Nation ^ 

In this paper, to make this statement more 
precise and systematic, I will apply the BRS gauge fixing procedure to the tensor model and 
numerically analyze the ghost kinetic term. Then I will compare the results of the numerical 
analysis with the corresponding BRS gauge fixing in the general relativity. 

This paper is organized as follows. In the following section, I will apply the BRS gauge 
fixing procedure to the tensor model. In Sec(3l I will discuss the corresponding BRS gauge 
fixing procedure in the general relativity. In SecJH I will numerically study the eigenvalues and 
eigenmodes of the ghost kinetic term in the tensor model at the Gaussian backgrounds with 
dimensions D = 1, 2, 3, and compare with the ghost kinetic term in the general relativity on 
the flat spaces in these dimensions. The final section is devoted to a summary and discussions. 



three-dimensional quantum gravity [201 |21j. See also [22j for more and the recent developments. 

* These modes of broken symmetries appeared as vanishing spectra of fluctuations in the previous works 

HEIZZ]. 

§The idea to consider local gauge symmetries to be non-linearly realized broken symmetries is rather old. 
For example, see [29l [30j [3l] . 
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2 BRS gauge fixing procedure in the tensor model 



2.1 Direct computation of the gauge volume 

Let me start with the direct computation of the gauge volume. 

The dynamical variable of the tensor model in this paper is given by a real rank-three 
tensor C a b c , which is totally symmetric, 

C a bc = Cbca = C ca0 = Cbac = Cacb = Ccba- (1) 

The index takes values 1, 2, . . . , N. There is also a nondynamical symmetric real tensor g ab , 
which is basically taken to be g ab = 5 ab . Therefore, by following the standard pairwise index 
contractions, the tensor model is invariant under the orthogonal group transformation O(N), 

C abc -> (MC) abc = M a a 'M h v M c c 'C a , h , d , (2) 

where M a a ' e O(N). This is the gauge symmetry of the tensor model. 

The O(N) symmetric metric in the space of the dynamical variable C is defined by0 

ds 2 c = dC ahc dC abc . (3) 

The inner product associated with the metric ([3]) between two rank-three totally symmetric 
tensors is defined by 

(A,B) = A abc B abc . (4) 
The infinitesimal SO(N) transformation of C is given by 

(rpi/~1\ rpi a! /~1 I rpi b 1 '/■"» i rpi d f~i fr\ 

y 1 < ~~')abc = 1 a ^a'bc-r 1 b ^ab'c + J- c ^abc'i \p) 

where T l a a ' (i — 1,2, ... , N(N — l)/2) are the real antisymmetric matrices forming the Lie- 
algebra so(N) in the vector representation. 

The volume measure in the space of C is defined from the metric ([3]). Dividing an in- 
finitesimal region into the gauge directions and the others, the infinitesimal volume dVc can 
be expressed as 

, ,rr± /Det((TC, TC)) 

dV c = dgdV+J V V ; - 6 

G y c V Det T,T ' v ; 



" There exists an ambiguity to add dC a t, b dC ac c to this metric. The addition will change some details of the 
analysis of both the tensor model and the continuum theory, but the mutual agreement should be obtained 
anyway. 
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where dg is the Haar measure of SO(N) and dV^r denotes the infinitesimal volume normal to 
the gauge directions. Here Det(- • • ) are the determinants of the matrices with components, 

(TC,TC) l3 = (TC,T*C), (7) 
(T, T)ij = /i Tr(W), (8) 

where Tr denotes the trace in the vector representation, and ho is a coefficient related to the 
normalization of the Haar measure. Since the integrand in ([6]) is invariant along the SO(N) 
gauge directions, the partial integration over the gauge directions is trivially performed as 



iVc = iVi y ° l{0{N) \rt {{ ^\ (9) 
SO(N) C n V Det « r . T » 

where n is the possible symmetry factor becoming larger than 1 if there exists a non-trivial 
M G SO(N) which satisfies C = MC. This factor n can practically be ignored, since the 
regions of such symmetric values of C have generally vanishing volumes in the space of C. 
Thus, ignoring all the factors independent of C, one finally obtains 



/ dV c = dV£ y/Det((TC, TC)). 

JsO(N) 



(10) 



'SO(N) 

2.2 BRS gauge fixing procedure in the tensor model 

I apply the general BRS gauge fixing scheme with the so-called B field presented in [32] to 
the SO(N) symmetry in the tensor model. The BRST transformation of C is given by 

(5 B c) abc = CjCrcOote, (ii) 

where are the ghost variables, which are assumed to be real. The BRST transformation of 
the ghost variables is given by 

$BC k = -PkCiCj, (12) 

where f 1 ^ is the structure constant of so(N), defined by [T l ,T^] = j n \T k . There are also the 
anti-ghost and the B- variables, the BRST transformations of which are given by 

5 B Ci = iBi. 

S B Bi = 0. (13) 

These c, and Bi are also assumed to be real. The nilpotency 5 B = can easily be shown by 
explicit computations. 
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The interest of the present paper is in the small fluctuations around certain backgrounds 
of C. Let me denote a background by C° and the fluctuations by A, 

C abc = C° abc + A abc . (14) 

Then the dynamical variable is shifted to A, and its BRST transformation is given by 

(S B A) abc = (5 B C) abc = Ci {TC) ahc = Ci {TC°) abc + Cl {TA) ahc . (15) 

The general scheme implies that the BRST exact action corresponding to the sum of the 
Faddeev-Popov and the gauge fixing terms can generally be given by 

S GF+FP = S B (ciF\A, c, c, B)) , (16) 

where F l are the (almost arbitrary) gauge-fixing functions with vanishing ghost number. A 
natural choice in the present case is 

F i = (T i C°,A), (17) 

since the gauge fixing conditions (F % = 0) only allow A to be normal to the gauge directions 
at the background C°. Computing ([TBI) with (fTTj) . Sgf+fp is explicitly given by 

Sgf + fp = iB^TC , A) - c^TC , T^C)c r (18) 

The path integral measure, which is just a usual integration in the present case, can be 
defined by 

[dA] Y[dBidcidci, (19) 

i 

where [dA] is the volume measure of A defined from the metric ds 2 A = dA abc dA abc , which is 
identical to the O(N) symmetric metric Q. Here a possible overall factor is not taken care 
of. From the O(N) invariance of the volume measure [dA], one can easily prove the BRST 
invariance of the integral, 

/ [dA] dBidadci 5 B {. . .) = 0, (20) 

i 

which guarantees the independence of physics from the choice of the gauge-fixing functions. 

2.3 Comparison between the direct and the BRS expressions 

In the following, let me compare the BRS result ([TB"]) . (fT9~j) with the direct computation ([TO]) . 
To do this, let me introduce a normalized orthogonal basis which divides the space about C° 
into the subspaces tangent {ff"} and normal {vf 1 } to the gauge directions, 

(rc% 0± > = o, 
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m / 


— fym 




.»?"> 


= o, 


k 011 







(21) 

In general, A can be expanded in terms of these vectors as 

A = a i v? +p l v'j ,± . (22) 

From the definition of the basis (l2~Tj) . the volume measure is [oL4] = P[ i da* d/? 4 . Putting 
( I22p into (1181) . and integrating over q, q, 5j and finally over a 1 , one obtains 



/ m n dB,^-^-^ = / n g^ff.^V " 8 ' '^'- (23) 

where an overall numerical constant is ignored, S(C) is the original unfixed action, and the 
matrices in the determinants are defined by 

{TC°,TC)ij = (rc°,T j c), 

(TCV 0|I >;; = (T*CV°">- (24) 

The result (|23|) does not look like (flOl) . but they are actually identical. To see this, let me 
introduce a similar normalized orthogonal basis {v\ }, {t>/~} around C = C° + A as fl2TT) . 

(T*C, U/" 1 ) = 0, 

(^^i 1 ) = 0, 

(wJ,«J) = ^, (25) 

Then the square of the determinants in (1231) can be computed as 

( SS^)| ) 2 = ((TC,TC ) ((TC ,f ")(f ",TC )) 1 (TC°,TC)^j 

= Det((TC,v°\)(v°\TC}) 

= Det ((TC,TC» [Det^'U 011 ))] 2 

= Det((TC,TC» [Det((^, W 0± ))] 2 , (26) 



where similar shorthand notations like (1241) are used to denote the matrices. In the above 
derivation, I have used the completeness of the bases {f "} , {t> "} in the subspaces tangent to 
the gauge directions, and 

[Det ({v«,,/'«))Y = Det ( ) = [Det ««\«"»]\ (27) 
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which can be shown from the identity, 

Det(D)Det(A - BD~ X G) = Det ( ^ D ) = Det ( A ) Det ( D ~ CA^B), (28) 

and the properties of the orthogonal normalized bases. From the definition (12"2"1) and that dV^ 
in fllOp is the infinitesimal volume normal to the gauge directions at C, one obtains 

dV£ = l[dp l | Det ((^,w 0± ))| . (29) 
i 

Thus (jTUj) and (1251) are actually identical, because of (1261) . (12"9~|) . 



3 BRS gauge fixing in the general relativity 

In this subsection, I will discuss the BRS gauge fixing procedure in the general relativity 
[331 [Ml EH [36] corresponding to that of the tensor model in the previous section. 

By rewriting the coordinate transformation of the metric tensor with the ghost fields, the 
BRST transformation of the metric tensor is given by 

&B9p» = V M cv + V v c^ (30) 

where V M is the covariant derivative, and c M is the ghost vector field. Then the nilpotency of 
the BRST transformation requires that the ghost field be transformed by 

(S B c^ = d>v v c? = c u d„c»). (31) 
The anti-ghost field and the S-field are introduced with the BRST transformations, 

5 B c^ = iB^, 

5bB^ = 0. (32) 

The nilpotency 5 2 B = can be checked by explicit computations^. All the fields above are 
assumed to be real. 

In the numerical analysis of the following section, I will take C° to be the Gaussian back- 
grounds [2"B"1 127] 127] , which correspond to the fuzzy flat spaces with arbitrary dimensions. 
Correspondingly, flat backgrounds are considered in the general relativity as 

Qiiv = <V + fyw> (33) 

"For example, 8 2 B c^, = can be shown from (5bT ^ v p = c° ' R lLa v p + V M V„c p and the Bianchi identities for 
the Riemann tensor. 
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where h^ v is the new dynamical field with = V M c^ + V u c^. 

In [27], it was argued that the metric ([3]) corresponds to the DeWitt supermetric 

ds] = j d D x^g-(g^g p ° + Ag^)5g^6g pa . (34) 

Thus the inner product associated to ( 134|) between two rank-two symmetric tensor fields is 
defined by 

(k, l) g = J d D x^(g^g^ + Ag^g™)k,J f „. (35) 

In the following, I want to obtain the BRS gauge fixing in the general relativity which is 
analogous to (|T6|) . (jTTj) . Since C° corresponds to the flat background, the analogy of T l C° 
are the infinitesimal local translations of the flat background. Therefore CjT l C° in the tensor 
model should correspond to the field d^c v + d v c^. The deviation A from the background in 
the tensor model corresponds to h^ u . Thus the action corresponding to ffT6l) . ffTTj) is obtained 
as 



^gf+fp — ^b{9ixC v + d v c v , h pcr ) g 

= 2 5 B (J d D x^(g^g^ + 4g*>gr)(d&)h 



(36) 



In the quadratic order of the fields around the flat background, ( 1361) becomes 

S g G % FP = 2 J d D x [i K + 4 d,B v hr) - 6 <9% d v c v - 4 d,c v . (37) 

The partial derivative of (1371) with respect to i?^ gives the gauge fixing condition as 

d^K + ±d V K„ = 0. (38) 

One can check that this gauge fixing condition is actually satisfied by all the metric fluctua- 
tion modes corresponding to the low-lying fluctuation modes in the tensor model which were 
reported previously in [271 EE] . 

Putting the form c M , c M = n^e ipx , the kinetic term of the ghost fields in (|37|) can be shown 
to have the spectra, 

20 p 2 for the longitudinal mode n M oc p M , , , 

8p 2 for the normal modes n^p^ = 0. 



Thus the longitudinal mode has no degeneracy, while the normal modes have the degeneracy 
D — 1 in D dimensions. 
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4 Numerical analysis of the ghost kinetic term in the 
tensor model 



In the papers [26J, [27J I2H], a class of tensor models possessing the classical solutions with 
Gaussian forms have been constructed and analyzed. In this section, I will take C° to be such 
Gaussian backgrounds, and numerically study the ghost kinetic term in the tensor model. 
In all the dimensional cases to be studied (D = 1,2,3), some massless trajectories of ghost 
modes, which are clearly separated from the other higher ghost modes, will be found, and they 
will be identified with the reparametrization ghosts in the general relativity. 



4.1 The coefficient matrix of the ghost kinetic term 

In the momentum basis, such a gaussian background C° has the form [261 [27J [28j [3] 

C pi,P2, P3 = ex P {- a (Pi +Pl+ pD) <Wp 2 +P3,o, (40) 

where a is a positive parameter, and each momentum is assumed to take integer vales bounded 
byL: 

D 

p=(p\p 2 ,...,p D ), 5>*) 2 <L 2 . (41) 

i=l 

Since the delta function in (1401) implies the momentum conservation, there remains the D- 
dimensional translational symmetry on this background. 

From (ITS]) , the ghost kinetic term is given by 

S gh = -c i (T i C°,TtC°)c j . (42) 

Because of the momentum conservation of the background (|40|) . it is most convenient to take 
a momentum basis for the generators T\ Namely, the indices of the generators are given by 
pairs of distinct momenta i — [pq] {jp ^ q), and the generators are antisymmetric matrices 
defined by 

(T^"^ S ) rs &p,r&q,s 3p,s$q,r- (43) 

Then, putting (T43!) into fT42|) . the coefficient matrix of the ghost kinetic term is given by 

[pigi],[p292] \ ' ' 

_ or r r r ,c \^ r r 

— OUp 1 J rP2 fiU qi J r q 2 fl / j >-v (?lir . jS >-' (?2j _ r , i _ s T U / j ^ 'p 1 ,52,r° 'qi,P2,-r 

r,s r 

-(pi ^ qx) - (P2 ^ q-i) + (pi ^ ?i,P2 ^ qi), (44) 
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where the summations of the momenta r, s are over the range (ETC]) , and the simplified notations 
for anti-symmetrization have been used. 

Because of the momentum conservation of the background, the matrix M? h , r , is di- 
vided into the block matrices with each value of the ghost momentum Pi + qx — — (P2 + 0.2)- 
Therefore, the analysis of M gh can be performed independently at each momentum sector. 

In the following subsections, I will study the spectra and the properties of the eigenmodes 
of M gh , for dimensions D = 1,2, 3, and compare with the continuum theory. 

The numerical facility was a Windows XP 64 workstation containing two Opteron 275 
processors and 8 GB memories. The C++ codeaE] were compiled by Intel C++ compiler 10.1 
with ACML 4.2 for Lapack/Blas routines. The output were analyzed in Mathematica 6.0. 

4.2 The D=l case 

In FigJU the eigenvalues of M 9h are plotted for two cases. A zero spectrum exists at p = 0, 
as expected from the fact that there remains a translational symmetry on the background. 
There clearly exist a series of spectra which form a massless trajectory and are clearly sep- 
arated from the other higher modes at low momenta. This series can be identified with the 
reparametrization ghost of the continuum theory, as explained below. In fact, the trajectory 
contains only one mode at each momentum value, which is consistent with the result fl39|) of 
the continuum theory. In the left figure with L — 15, a — 0.5/ L 2 , the trajectory looks to 
have a linear momentum dependence near the origin, which contradicts (I3"§|) . But, as can be 
seen in the right figure with L = 1500, a = 3/L 2 , the trajectory tends to become reasonably 
smooth near the origin in the cases with larger L and a, which is consistent with ( |39|) . This 
is in agreement with the natural expectation that the continuum theory will be obtained only 
in large L and at low momenta. 

4.3 The D=2 case 

The result (|39|) of the continuum theory implies that there should exist two massless trajec- 
tories of spectra with no degeneracy, and that the ratios of the spectra in the two trajectories 
should be given by 22 =2.5. In fact, in the left figure of Figj2j one can find that there exist two 
massless trajectories linked to the two zero spectra at p — 0, which come from the unbroken 
translational symmetries. The numerical data also show that each trajectory has only one 
mode at each momentum value. The right figure shows that the ratios of the two trajectories 
at each momentum are actually in good agreement with 22. 

**The codes are downloadable from http://www2.yukawa.kyoto-u.ac.jp/"'sasakura/codes/ghostcpp.zip 
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Figure 1: The eigenvalue plots for D — 1. The horizontal axis is the momentum of the ghosts. 
The left figure shows the whole spectra for L — 15, a — 0.5/ L 2 . The right figure shows the 
low part of the spectra for L = 1500, a = 3/L 2 . The fitting line is 8.1 x 10~ 4 p 2 — 2.1 x 10~ 7 p 4 . 



' .s 



R,ltK> 

3.0 r 



Figure 2: The left figure shows the low part of the spectra for D = 2, L = 100, a = 2/L 2 . The 
right figure shows the ratios of the two trajectories. The horizontal axis is the momentum size 
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Figure 3: The contour plots of 5K q _ q+p for the modes contained in the two massless trajectories 
for D = 2, L = 20, a = 2/L 2 . The horizontal axes are (q , q 2 ) = q. The left and right figures 
are shown for the modes in the lower and upper trajectories, respectively, at the momentum 
p=(l,0). 

It will also be a good check to see whether the mode profiles are consistent with the 
continuum theory. To do this, I follow the strategy in the previous works [27J [28] . Let me 
define a tensor, 

K ab = C acd C b cd . (45) 
Under small fluctuations around C°, this tensor fluctuates as 

5K ab = 5C acd C\ cd + C° acd 6C b cd . (46) 

The present interest is in the fluctuations in the gauge directions T l C°. For the gauge direction 
determined by an eigenvector v of M 9h , 5K is given by 

6K v ab = {viTCX^C ^ + C\ cd { Vi rC\ cd . (47) 

In Figj3l 5K is plotted for the eigenmodes contained in the two trajectories. 

On the other hand, the correspondence between the tensor model and the general relativity 
implies $27\ EE] 

5K PlP2 = h lu ,(p 1 +p 2 ) (pi-P2T(pi -p 2 yexp ~P2) 2 ^J ■ (48) 

Substituting the gauge transformation h fiv (p) = n^pv + n u p u into this expression, one obtains 
the two contour plots in FigJUfor the normal (p^n^ 1 = 0) and the longitudinal (p^ oc n M ) modes, 
respectively. These figures are in full agreement with Figj3] in view of the mode assignment 

Q32). 
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Figure 4: The contour plots of 5K q - q+p for p = (1, 0) expected from the continuum theory. 
The left and right figures are for the normal = 0) and the longitudinal oc n^) modes, 

respectively. 
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Figure 5: The left figure shows the low part of the spectra for D 
The right figure shows the ratio of the two trajectories. 



3, L = 15, a = 1.5/ L 2 



4.4 The D = 3 case 

The result (1391) of the continuum theory implies that the lower trajectory should contain two 
modes at each momentum. In the left figure of FigJSJ the low part of the spectra for D = 3, 
L = 15, a = 1.5/ L 2 is shown. There exists two massless trajectories, and in fact the numerical 
data show that the lower trajectory contains two modes at each momentum value. In the right 
figure, the ratios of the two trajectories are shown, which are consistent with ^. 
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5 Summary and discussions 



In this paper, I have applied the BRS gauge fixing procedure to the tensor model, and have 
numerically analyzed the ghost kinetic term at the Gaussian backgrounds, which correspond to 
the fuzzy flat spaces with arbitrary dimensions. Then it has been found that there exist some 
massless trajectories of the ghost modes, which are clearly separated from the other higher 
ghost modes. By examining the properties of the modes in these massless trajectories, it has 
been shown that these modes can be identified with the reparametrization ghosts in the BRS 
gauge fixing of the general relativity. This means physically that the local gauge symmetry 
(the local translation symmetry) is emergent around these backgrounds in the tensor model. 

Combined with the results of the previous works, this paper has shown that the low-lying 
fluctuations around the Gaussian backgrounds in the tensor model correctly generates the 
general relativity, including its gauge symmetry, on the flat spaces in general dimensions. 
However, this has only been shown in the quadratic order of the fluctuations around the 
backgrounds, but not for any of the higher non-linear terms. On the other hand, the general 
relativity (possibly with modification of the action) is the unique interacting field theory of 
the rank-two symmetric tensor field with the gauge symmetry. Therefore there exists a good 
chance for the tensor model to correctly generate also the non-linear terms. This should be 
studied in future works. 

The above agreement between the tensor model and the general relativity including the 
gauge symmetry also suggests that the quantization of the general relativity can be realized 
by that of the tensor model, and thus a kind of quantum gravity can be defined by the 
tensor model. There exist a lot of questions to be addressed by quantum gravity, the most 
phenomenologically interesting of which would be the cosmological constant problem [38]. In 
the conventional approaches, one needs an extreme fine-tuning of the cosmological constant to 
stabilize a flat space from quantum corrections. Therefore it would be interesting to study how 
the Gaussian backgrounds, which represent fuzzy flat spaces, respond to quantum corrections 
in the tensor model. 

So far, the correspondence between the tensor model and the general relativity has been 
shown only for a limited class of tensor models, which have the Gaussian solutions. It is 
also an interesting future problem to investigate whether such correspondence holds in more 
general tensor models. 
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